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Abstract. The two dimensional crossover from independent particle towards collective motion is studied 
using 2 spinless fermions interacting via a U/r Coulomb repulsion in a L x L square lattice with periodic 
boundary conditions and nearest neighbor hopping t. Three regimes characterize the ground state when 
U /t increases. Firstly, when the fluctuation Ar of the spacing r between the two particles is larger than 
the lattice spacing a, there is a scaling length Lq = VSn 2 (t/U) such that the relative fluctuation Ar/{r) 
is a universal function of the dimensionless ratio L/Lo, up to finite size corrections of order L~ 2 . L < Lq 
and L > Lo are respectively the limits of the free particle Fermi motion and of the correlated motion of 
a Wigner molecule. Secondly, when U/t exceeds a threshold U*(L)/t, Ar becomes smaller than a, giving 
rise to a correlated lattice regime where the previous scaling breaks down and analytical expansions in 
powers of t/U become valid. A weak random potential reduces the scaling length and favors the correlated 
motion. 
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1 Introduction 

An important issue in quantum many body theory is to 
know how one goes from independent particle motion to- 
wards collective motion when one decreases the density 
n s of electrons repelling each other via a U/r Coulomb 
repulsion. One can trap the electrons by a positive back- 
ground of charges (jellium model) unless one uses periodic 
boundary conditions (2d torus geometry). Many studies 
|jl|,||J|^,^,|U^] assume a parabolic trap, its rotational in- 
variance allowing to decouple the motion of the center of 
mass from the relative motions. A parabolic confinement 
has the merit to be realizable using semiconductor het- 
erostructures HH] for electrons or electromagnetic fields 
for cold ions. It has however the disadvantage to yield a 
non uniform charge density, so that the formation of the 
electron solid when the confinement becomes weak results 
from a complicated interplay between edge and bulk or- 
derings @||. This leads us to study a system without edge, 
having the geometry of a two dimensional torus, where the 
density is uniform. For a continuous torus, the size of the 
corresponding Hilbert space is infinite, making the use of 
a truncated basis unavoidable for a numerical study. A 
way to avoid such a truncation is to take a lattice model 
(tight binding approximation). Considering the the sim- 
plest limit: the ground state (GS) of two spinless fermions 
in a square L x L lattice with periodic boundary condi- 
tions (BCs), we give an answer to four questions. Is there a 
simple one parameter scaling theory of the Fermi- Wigner 



crossover? How important are the finite size corrections 
to this scaling theory? Does the lattice play an impor- 
tant role? How those answers are modified when a weak 
random potential is included? 

The answers can be summarized as follows. In the 
strong coupling limit, the Coulomb repulsion pushes the 
spacing |r| between the two particle to take its maximum 
value L/y/2 and its fluctuations Ar become of the or- 
der of the lattice spacing a. This gives rise to a strongly 
correlated lattice regime where analytical expansions in 
powers of t/U are sufficient to describe the system. For 
weaker couplings, the relative motion becomes more and 
more delocalized and the lattice effects become irrelevant. 
This gives rise to a universal regime characterized by a 
scaling length L = \/8n 2 (t/U). This length corresponds 
to the system size for which the extra energy to add a sec- 
ond particle in the system is the same in the Wigner limit 
(t = 0) than in the Fermi limit (U = 0). L < L is the 
Fermi limit, L > Lq is the Wigner limit and the Fermi- 
Wigner crossover takes place when L = Lq. In this uni- 
versal regime, the relative fluctuations Ar/ (r) = F(L/Lg) 
up to finite size corrections of order L~ 2 . When the site 
potentials have weak random fluctuations of order W, 
there is still a regime where one parameter scaling re- 
mains valid, Ar/(r) being characterized by a function 
F\y(L/Lw) which depends on W. Both Ar/(r) and the 
scaling length Lyy become smaller, in agreement with the 
general idea that a weak disorder favors the correlated 
motion. 
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2 Lattice model 

We consider two electrons with symmetric spin wave func- 
tions and antisymmetric orbital wave functions (spinless 
fcrmions), free to move in a L x L lattice with periodic 
BCs, and interacting via a U/\r\ repulsion. The Hamilto- 
nian reads 



fl = -t V (cjc.,- + h.c) + y^Vj 

<i,j> i 



u 



E 



2 ir; 



(1) 



where i, j label the lattice sites, < i,j > means i near- 
est neighbor to j, cj, Ci are the creation, annihilation op- 
erators of a spinless fermion at the site i; rn = cjcf is 
the occupation number at the site labelled by the vector 



i = (i Xl i y ). The vector 



j is defined as the short- 



est vector going from the site i to the site j in a square 
lattice with periodic BCs. This means that r x < L/2 and 



f y < L/2 if L is even, L 



1 if L is odd. Hence, 



the pairwise interaction U/\v\ exhibits a singularity on 
the lines r x — L/2 and r y = L/2 and another one at 
their crossing point r x = r y = L/2. t = h 2 /(2ma 2 ) is 
the hopping term, a the lattice spacing, Vi the site po- 
tentials which arc randomly distributed in the interval 
[-W/2,W/2] and U = e 2 /(ea) the Coulomb interaction 
between two fermions separated by a in a medium of di- 
electric constant e. When there is no disorder (W = 0), 
k = (k x ,k y ) being the one particle momentum, it is more 
convenient to write H using the Fourier transforms of the 
creation and annihilation operators. One has the relations 



T.T, 



a k e 



ikj 



and 



(2) 



(3) 



for N = 2 and M = L 2 . When L is odd, this gives M 
blocks of dimension (M — l)/2. When L is even, we have 
two different dimensions of the blocks. 



N H = M bl M 1 + M b2 M 2 



(8) 



where Mi , M 2 are the number of blocks with dimensions 
M bl , M b2 . Here, 

Mi = M ^~? H , M 2 = N «- M ^ (9) 



M b2 - M bl ' 



M h2 - M bl 



and 



M n M 
M bl = — -2, M b2 = —. 



(10) 



3 The Fermi limit 



When U = W = and with periodic BCs, the states are 
Nh plane wave Slater determinants 



«|0) = |k 1 k 2 ). 
k x = (2-n/V)n x and k y = (2ir/L)n y with 



0,±1,...,±±^± for L odd 

;/ • "" ~ ^ 0,±1,...,±(|-1)4 for L even 

The two particle wave functions 



V>(ri,r 2 ) = -j= 



l.gjki-ri _^ e iki-r 2 



J_ e ik 2 ri J_ e ik 2 r 2 



become 



i'(ri.n) = ( y iKR ) (i^sink-r 



(11) 



(12) 



(13) 



(14) 



which yield 



after the change of coordinates: 



ff = E a k fl k £ ( k )+ E a k 2 +q a k 2 4 1 -q a k 1 ^(q) (4) 
k q,ki,k 2 



where 
and 



e(k) = — 2t (cos k x + cos k y ) 



V{ci) 



E 



iqr 



2L 2 ^ r 

r#(0,0) 



(5) 
(6) 



Without disorder, the total momentum K = ki +k 2 is 
conserved. The two particles states of different momenta 
are not coupled, H becoming block diagonal, each block 
corresponding to a given total momentum K. 

The dimension of the two particle Hilbert space is 
given by 



N 



Ml 



M — 1 



H 



N\(M — N)\ 



-M 



(7) 



K = ki+k 2 , 



R = -(n+r 2 ) 



k=-(k 2 -kx) 



r 2 - ri 



(15) 
(16) 



Without disorder, the motion of the center of mass and 
the relative motion are separable: 



V(ri,r 2 ) = 0(R) X (r). 



(17) 



(f> (R) is a plane wave of total momentum K which de- 
scribes the propagation of the center of mass on the 2d 
torus while \ ( r ) describes the relative motion. The mo- 
ments (|r| m ) of the inter-particle spacing r are given by: 



(|rD 



|r| m p(r)dr 



(18) 



where p(r) = |x( r )| 2 is the inter-particle spacing distribu- 
tion. 
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The one particle energies e(k) can be ordered by in- 
creasing values. The one particle GS energy £o is not de- 
generate, while the three next excitations £1, £2, £3, are 
four-fold degenerate; the fourth excitation £4 is ten-fold 
degenerate; etc. 



£0 
£1 

£2 

£3 

e 4 



-At 

-2t I 1 + cos — 

-4f cos — 
L 

4-7T 

-2i ( 1 + cos — 



, 2n 47T 
-2f ( cos — -(- cos — 



(19) 



p(rx,ry) 



0.0004 




Fig. 1. Inter-particle spacing distribution p(r x , r y ) for L = 60, 
[7 = and K = (0,2tt/L). 



The two particle GS consists of one particle of energy 
£0 and of a second particle of energy e\. Because £1 has a 
four-fold degeneracy, the two particle GS energy 



E = £0 + £1 = -4i - 2* 1 



2tt 



cos ■ 



(20) 



is also four-fold degenerate. 

Hereafter, we study the two particle GS of momentum 
K = (0,2tt/Z). For ki = (0,0) and k 2 = (0,^), the GS 
wave function is given by 



( Mr,.r_.) = ( ^e l ^ R y 



which yields 



.a/2 . 7rr.„ 
-%— — sin — - 
Li ij 



P(r x ,r y ) = — sin — 
4 E (^+^) m/2 sin 2 ^. 



(21) 

(22) 
(23) 



The corresponding inter-particle spacing distribution 
p(r x , r v ) is shown in Fig. Q for L = 60. 

When one turns on the interaction U, the GS wave 
function of total momentum K = k; + kj can be written 

as 

\MU))=J2 a ^ ( f/ ) I k > k i ) ■ ( 24 ) 

i<j 

The coefficients a^k^U) — (kjkj IV'o(C^)) can be numeri- 
cally obtained if L is not too large, and the inter-particle 
spacing distribution p(r x ,r y ) given by 

Z J j<j,m<n tt kik3 a k m k n o k . +kj |_e e J 

is shown in Fig. | for L = 60, = 5 and K = (0, 2%/L). 
If one particle is located at the site (0,0), one can see 
how the interaction localizes the second one near the site 
r = (i/2,i/2). 



p(rx,ry) 
0.004 



0.002 




Fig. 2. Inter-particle spacing distribution p(r x ,r y ) for L = 60, 
U/t = 5 and K = (0,2tt/L). 

4 The correlated lattice limit 

Before studying in more details the interaction induced lo- 
calization of the inter-particle spacing r, let us consider the 
other limit t/U — > where the two particle wave functions 
\tp) are more conveniently written using the operators cj. 
In this limit, the kinetic part 



H 



K 



= -*£ * 



(25) 



<i]> 



of .ff is a small perturbation compared to the Coulomb 
part and the levels can be expanded in powers of t/U. 
When t/U — ► 0, at the order zero of a t/U expansion, 
the Coulomb energies Eq and E\ of the GS wave func- 
tion (K)) and of the two degenerate first excitations 
|?/4°^(K)) and ^^(K)) are given by 



and 



E l 



E 



U_ 



U 
do 



V2U 



u 



v(t) 2 +(f-lr 



(26) 
(27) 



4 



M. Martinez and J.-L. Pichard: Independent particle and collective motion 



if L is even. When L is odd, the lengths do and d-y must 
be differently defined: 



d = 



L-l 

~7T 



and 



di 



Y (I. I i x 2 



(28) 



(29) 



Hereafter, we assume that L is even and we con- 
sider the states of total momentum K = (0,2%/L). a = 
(L/2, L/2). bi = (L/2 + 1, L/2) and b 2 = (L/2, L/2 + 1) 
defining the locations of three lattice sites which are as 
far as possible from the site = (0, 0), the zero order GS 
wave function 



l4 0) (K)) = -^^exp^cH 4 



JO) (30) 



is directly coupled at the order t/U to the two wave func- 
tions of energy E\ 

|< ) (K)) = i^exp jK -J c ] c ] +bi |0) (31) 



and 



by 



and 



|4 0) (K)) = i^exp^ c ] c ] +b2 |0). (32) 



(4 0) (K)\H K \4 0) (K)) = -2V2t 



(33) 



(4 0) (K)\H K \4°\K)) = -V2 1 + (34) 



2tt 



L 



respectively. At the order t/U, the GS wave function is 
given by 

^^^|V,r(K)) + |^ 1) (K)) (35) 



where C 2 is a normalization constant and 

± ^ 0) < K)\H K \4 0) < 
E — Ei 



(36) 



a=l 



One gets 



\S4 X) (K))=A\^>(K))+B\^(K)) (37) 



s(0) 



(0), 



where 



A = 2^2 



d$di ( t 



do — di \U 



(38) 



and 



(40) 



-1.0 



-1.5 - 



-2.0 - 



1 

U/t 



Fie. 3. Relative fluctuation function of U for L = 20 

and K = (0,2-k/L). The log-log plot shows the three regimes 
characterizing the ground state. The three dashed lines corre- 
spond to u r — 0.25, u r oc (i/i7) ' 31 and the t/U perturbative 
behavior given by Eq. (|42|). 



for L even and K = (0, 2tt/L). This yields for the fluctu- 
ation At = y/ (r 2 ) — (?*) 2 of the inter-particle spacing: 



Ar=^LJ[(L-iy + l}U + cos^) (41] 



and for its relative fluctuation u r = Ar/(r) 

u r = A J 2[{L - l) 2 + 1] (3 + cos ^ ) ( 42 i 



2tt 



This correlated lattice behavior is shown in Fig. || for 
L = 20 when L//t > 1000. 



5 The three regimes of coupling 

The behavior of the relative fluctuation u r = Ar/ < r > 
is shown in Fig. | for L = 20 and K = (0,2tt/L). 
The crossover from independent particle towards strongly 
correlated motion in a finite size lattice exhibits three 
regimes: 

— A weak U Fermi regime where u r is large and almost 
independent of U. The fluctuation Ar is of the order of 
the expectation value (r), r being broadly distributed. 

— An intermediate correlated Wigner regime where u r 
exhibits a weak (t/U) a decay, where a w 0.31. Ar 
becomes small compared to (r), the distribution of r 
becoming narrower. 

— A large U strongly correlated lattice regime where u r 
exhibits the strong t/U decay predicted by the pertur- 
bative expansion (Eq. (fl2"|)), the distribution of r being 
extremely narrow. 
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Fig. 4. Relative fluctuation u r — g(L) where g(L) — 0.9/(L + 
l) 2 as a function of L/Lq w UL/(\/ r 8iv 2 t) for even L. The log- 
log plot shows the universal scaling function F(L/Lq) for U < 
U*(L), and the non universal lattice regime for U > U*(L). 
The filled circles approximately give the thresholds U*(L). 

The value U*(L)/t characterizing the crossover be- 
tween the intermediate Wigner regime and the strongly 
correlated lattice regime is consistent with the condition 
Ar < 1. When U ~ U*(L), the lattice strongly reduces 
the degrees of freedom of the relative motion, and lattice 
expansions in powers of t/U become valid. We study in 
the next section the Fermi- Wigner crossover which occurs 
before this strongly correlated lattice regime. 



6 Scaling theory of the Fermi-Wigner 
crossover 

For U < U*(L), we define the scale Lq associated to the 
interaction strength U/t such that a dimensionless observ- 
able as u r becomes a universal function of the dimension- 
less ratio Lj Lq, up to certain finite size corrections which 
we will estimate. The argument for defining the scaling 
length Lq can be presented as follows. Let us consider the 
system at U = where the first particle occupies the state 
£o = —At. Due to Pauli principle, adding a second particle 
requires an extra kinetic energy: 



AE(U = 0) = ex - e = 2t(l - cos(2tt/£)) 



(43) 



In the other limit t — ► 0, adding a second particle requires 
an extra Coulomb energy: 



AE(t = 0) = 

do 



(44) 



The length Lo characterizes the scale at which 
AE(U = 0) coming from Pauli exclusion principle is equal 
to AE(t = 0) coming from Coulomb repulsion. When L 
is large enough, Lq can be approximated by: 



(45) 



3000 



2000 - 



1000 




Fig. 5. Circles: Interaction thresholds U*(L) giving the end of 
the universal scaling regime and the beggining of the strongly 
correlated lattice regime, extracted from Fig. |^. The dashed 
line corresponds to the condition Ar(U) = 2/3 where Ar is 
given by Eq. (fll|). 



but when L is small, Lq should be more precisely deter- 
mined by solving the equation AE(U = 0) = AE(t = 0), 
and without neglecting the even-odd effects in the defi- 
nition of do- For scales L < Lo, the GS wave function 
has mainly to minimize the kinetic energy while the mini- 
mization of the Coulomb energy becomes more important 
for the scales L > Lq. The Fermi-Wigner crossover is ex- 
pected at L/Lq = 1. 

One assumes for the dimensionless ratio u r the usual 
finite size scaling ansatz [ [To[ : 



u r (L, U, t) 



F(±-)+g(L), 



(46) 



the finite size correction g(L) 
be easily evaluated when U/t 
becoming 



> as L — > oo. g(L) can 
(Lq — > oo), the ansatz 



.{L,U = 0,t) = F(0)+g(L) 



(47) 



This is a problem without interaction where one can use 
Eq. (p|), for eventually obtaining F(0) = 0.2543 and 



J(L) 



0.9/ (L + lf fori even 
-0.4/ (L - l) 2 for L odd 



(48) 



The finite size correction g(L) has a sign which depends 
on the parity of L and disappears as the inverse of the 
total number of lattice sites. 

Fig. |] gives F(L/L ) — u r — g(L) for even values of L. 
One can see the universal scaling behavior with a Fermi 
Wigner crossover at L = Lq up to a value U* (L) denoted 
by the filled circles where scaling breaks down and the 
correlated lattice regime begins. The behavior of U*(L) is 
shown in Fig. S, where one can see that the values U* (L) 
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extracted from Fig. ^ are precisely given by the condition 
Ar S3 2/3, Ar being given by Eq. @. Below U*(L), 
u r — g(L) is a one parameter function of L/Lq oc LU/t, 
and the lattice effects are irrelevant. The lattice regime 
occurs only above a large interaction U*/t oc 1? when L 
is large. 

For N spinless fermions in a L x Z, square lattice, the 
Coulomb energy to Fermi energy ratio becomes pM: 



£7 1 



2f 



(49) 



for a density i/ = N/L 2 . As it must be, the dimension- 
less ratios L/Lq and r s are proportional, the advantage of 
using L /Lo instead of r s being that the crossover occurs 
when L/Lq = 1, and not at r s = 1. In the strongly corre- 
lated lattice regime, u r becomes a function of Lt/U , and 
not of r s oc LU jt. 

The Fermi- Wigner universal crossover is illustrated 
with more details in Fig. |^ where we have only taken 
even values of L. One can see that the finite size correc- 
tion g(L) = 0.9/(L + l) 2 allows us to plot all the values 
u r — g(L) calculated for U < U*(L) onto the universal 
scaling curve F(L/L ) when L = 6, 10, . . . , 60. Fig. [7| gives 
the same universal curve F(L/Lq) for the odd values of 
L, once the finite size correction g(L) = — 0A/(L + l) 2 
have been subtracted to u r . The universal function F(x) 
is close to its value F(0) = 0.2543 when x < 1 and behaves 
as x~ a with a S3 0.31 when x > 1. 

The value a S3 0.31 depends on the exact form of the 
pairwise repulsion U/\r\ around r = (L/2, L/2). Due to 
our definition of the inter-particle spacing r, U(r) is non 
analytic around r = (L/2, L/2), making the study a little 
bit involved. If one modifies the pairwise repulsion such 
that U(r) becomes analytic around (L/2, L/2), taking for 
instance 



U(r) 



U 



sin 2 (¥)+sin 2m 



(50) 



instead of U/\r\, one can show |l2| that the equation for 
the relative motion becomes identical to the single particle 
motion in a 2d harmonic potential when U is large enough, 
to eventually obtain a slightly different exponent a = 1/4 
in the continuous limit (L — > oo). 




Fig. 6. Fermi- Wigner crossover for the relative fluctuation u r , 
described by the universal scaling function F(L/Lq) = u r — g(l) 
in a log-log plot. Values calculated from even values of L with 
the finite size correction g(L) — 0.9/(L + l) 2 . The solid lines 
correspond to F(x) — 0.2543 and F(x) oc a; -0,31 respectively. 



fa. 




-0.9 - 



Fig. 7. Fermi- Wigner crossover for the relative fluctuation u r , 
described by the universal scaling function F(L/Lq) — u r — 
g(L) in a log-log plot. Values calculated for odd values of L 
with the finite size correction g(L) = —0A/(L — l) 2 . 



7 Effects of a weak random substrate 

We now study the effect of random site potentials Vi upon 
the Fermi- Wigner crossover. We restrict our study to the 
limit of weak disorder, the values of W being large enough 
to yield one particle diffusive dynamics in the non interact- 
ing system, but too small to yield one particle Anderson 
localization. In the presence of disorder, the momenta K 
are not conserved, and we need to diagonalize the total 
Hamiltonian H given by Eq. (|l|) in the site basis c|cj |0). 
Moreover, the relative fluctuations u r are now random 
variables, and the ensemble average values (u r ) will be 
considered. 



Figure g shows for L = 14 the difference between 
the u r obtained for W = and the (u r ) calculated for 
1 < W < 4. One can see that the disorder reduces the rel- 
ative fluctuations (u r ), shifts the Fermi- Wigner crossover 
to weaker interactions, and reduces the exponent a char- 
acteristic of the U~ a decay for intermediate interactions. 

For a weak value of W and small enough values of 
U/t to avoid the correlated lattice regime, we assume the 
generalized scaling ansatz : 

(u r (L, U, t, W)) = F w ( f ) + g w {L). (51) 
L w (U,t, W ) 
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Fig. 8. Ensemble average relative fluctuation (u r ) as a func- 
tion of U/t for L = 14 and W = 0, 1, 2, 3, 4 in a log-log plot. 



-0.6 




-1.4 1 1 1 1 1 1 1 1 1 1 1 

-3-2-10 1 2 

Fig. 9. Scaling function Gw{L/Lw) + C for W = 2 and 
6 < L < 14. An arbitrary constant C = 0.2 has been taken for 
the log- log plot. 

Since 

ML, U = 0,t, W)) = F w (0) + g w (L), (52) 

which implies that 

G w = {u r (L, U, t, W)) - (u r (L, U = 0,t, W)) (53) 

must only depend on the ratio L/L\y- The scaling length 
Lw is defined in such a way that the Fermi-Wigner 
crossover remains at L/Lw = 1. As shown in Fig. all 
the data calculated for L = 8, 10, 12, 14 can be mapped 
onto the same scaling curve Gw{L/ Lw) = Fw{L/Lw) — 
-FV(O). We just show the curve Gw=2(L/Lw), the curves 
obtained for W — 1,2,3,4 being also consistent with the 
ansatz ( |5l| ) when U and L are varied in the same range. 
G w (x) w for x < 1 while G w (x) w for i > 1, 



0.32 




0.6 1 ' 1 1 1 1 1 1 1 

12 3 4 

W 



Fig. 10. As a function of W, exponent a(W) characterizing 
the the scaling function Gw{x) oc x~ a( - w ^ when x > 1 (upper 
figure) and reduction factor Lw / Lo for the scaling length Lw 
(lower figure). 



the exponent a(W) being given in upper Fig. [T^. To obtain 
that the Fermi-Wigner crossover remains at L/Lw = 1, 
one needs to take the scaling lengths L\y, the reduction 
factors L\y/Lo being given in the lower Fig. [10[ 

The reduction of the crossover scale Lw by the ran- 
dom potentials can be qualitatively explained if we revisit 
the argument which has allowed us to define Lq with- 
out disorder: L = Lw when the energies AEw which are 
necessary for adding a second particle in the one particle 
system are the same in the two limits U = and t = 0. 
On one hand, when U = 0, the first one particle excitation 
has a fourfold degeneracy which is removed by disorder. 
This gives a reduction of the energy cost (AEw{U = 0)) 
necessary in the Fermi limit. On the other hand, the en- 
ergy cost (AEw{t = 0)) necessary in the Coulomb limit 
remains almost unchanged if W is weak. This gives a shift 
of the Fermi-Wigner crossover to a lower scale Lw < Lq 
in the limit of a weak disorder. In Fig. |ll|, we have plotted 
as a function of L for W — and W = 2 the disorder 
averaged energies (AE W (U — 0)) and (AE w {t = 0)). 
One can see indeed that (AEw{U = 0)) is decreased by 
the disorder, while (AEw(t = 0)) remains essentially un- 
changed. The curves cross at L = Lo w 8.9 when W = 
and L — Lw ~ 7.1 when W — 2. The obtained ratio 
Lw=2/Lq w 0.8 is indeed smaller than 1, though some- 
what below the value Lw=2/Lq ~ 0.9 given in Fig. [n]. 

This quantitative disagreement may be due the distri- 
butions of the considered random variables. We have con- 
sidered only the ensemble average behaviors, which may 
not give the typical scaling behaviors, both for the rel- 
ative fluctuations and for the energies. When U = 0, it 
is indeed well known that the level spacing distributions 
are not normally distributed, but depend on the nature of 
the one particle dynamics. In the bulk of the spectrum, 
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10. 



11. 

12. 
13. 



5. 
6. 

and P. Hawrylak 
7. A.V. Filinov, M. 



J. Kyriadidis, M. Pioro-Ladriere, M. C ioga, A.S. Sachrajda 



:ond-mat/01 11543 



Bonitz and Yu.E. Lozovik, Phys. Rev. 
Lett. B 86, 3851 (2001). 

L. Jacak, P. Hawrylak and A. Wojs, Quantum dots, 
(Springer- Verlag, Berlin Heidelberg, 1998). 
I.V. Schweigert, V.A. Schweigert and F.M. Peeters, Phys. 
Rev. Lett. 84, 4381 (2000). 

M.E. Fisher and M.N. Barber, Phys. Rev. Lett. 28, 1516 

(1972); J.-L. Pichard and G. Sarma, J. Phys. C: Solid State 

Physics 14, L127 and L617 (1981) and refs. therein. 

G. Benenti, X. Waintal, and J.-L. Pichard, Phys. Rev. Lett. 

83, 1826 (1999). 

Z. A. Nemeth, unpublished. 

S.T. Chui and B. Tanatar, Phys. Rev. Lett. 74, 458 (1995). 



In conclusion, let us underline that an Hamiltonian 
with three parameters, U, t and W, has many different 
limits. This section is restricted to the study of the com- 
petition between the interaction U and the kinetic energy 
t when the disorder W remains a weak perturbation, in 
the limit L < L±, L\ denoting the one particle localiza- 
tion length. In this limit, we have shown that (u r ) is still 
given by a function F\y of the dimensionless ratio L/ Lw 
and that the scaling length L\y decreases as a function 
of W . This is in agreement with the general idea that a 
random substrate favors charge crystallization [[ll|]l3| , the 
correlated motion persisting to larger densities (weaker 
values of r s ). But for larger W, the issue is to describe the 
crossover from a Fermi glass of Anderson localized states 
towards a Wigner glass. In this strongly disordered limit, 
the u r will have large sample to sample fluctuations, and 
one cannot rule out that the average (u r ) will finish to 
be meaningless or that the one parameter scaling func- 
tion Fw{L/Lw) could become a more complicated two 
parameter scaling function of L/Lq and L/L\. The study 
of those glassy behaviors is left for another study. 



